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Abstract. An approach based on the spectral and Lie - algebraic techniques for constructing 
vertex operator representation for solutions to a Riemann type hydrodynamical hierarchy is 
devised. A functional representation generating an infinite hierarchy of dispersive Lax type 
integrable flows is obtained. 



1. Introduction 

Nonlinear hydrodynamic equations are of constant interest still from classical works by B. Rie- 
mann, who had extensively studied them in general three-dimensional case, having paid special 
attention to their one-dimensional spatial reduction, for which he devised the generalized method 
of characteristics and Riemann invariants. These methods appeared to be very effective [T] in 
investigating many types of nonlinear spatially one-dimensional systems of hydrodynamical type 
and, in particular, the characteristics method in the form of a "reciprocal" transformation of vari- 
ables has been used recently in studying a so called Gurevich-Zybin system [51 [3] in [5] and a 
Whitham type system in [S1[S]. Moreover, this method was further effectively applied to studying 
solutions to a generalized (owing to D. Holm and M. Pavlov) Riemann type hydrodynamical 
system 

(1.1) D^u^O, Dt:=d/dt + ud/dx, 

where N G Z+ and u G C°°(M'^; M) is a smooth function. Making use of novel methods, devised in 
[24117] and based both on the spectral theory [9l [TBI [191 [18] and the differential algebra techniques, 
the Lax type representations for the cases = 1,4 were constructed in explicit form. 

In this work we are interested in constructing a so called vertex operator representation |131 1141 
[?T1 ?] for solutions to the Gurevich-Zybin hydrodynamical hierarchy (|l.ip at = 2 : 

DfU — Uf + UUx = V, 

Dtv = vt+ uvx 0, 



(1.2) 

making use an approach recently devised in [221 [23] for the case of the classical AKNS hierarchy of 
integrable flows, and which can be easily generalized for treating the problem for arbitrary integers 

N e Z+. 

2. A VERTEX OPERATOR ANALYSIS 



(2.1 



We begin with a Lax type linear spectral problem [3 [4] for the equation (|1.1[) at = 2 

DfU — Uf + UUx = V, 

DtV = vt + uvx = 0, 



defined on the space of smooth real- valued 27r-periodic functions (w, w)^ € M C C°°(R/27rZ; R^) 
(2.2) df/dx= e[u,v;X]f, ^l^'^'>^]--=('xy2^ Xu7/2 
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where, by definition, v :— DfU, f E Loo(IR/27rZ; C^) and A G C is a spectral parame- 
ter. Assume that a vector function {u,v)^ e M depends parametrically on the infinite set 
t :— {ti,t2,t3, . . .} E ]R^+ in such a way that the generalized Floquet spectrum [9l [151 HI] 
a{£) := {A G C : sup3.gR A)||oo < 00} of the linear problem ()2.2|) persists in being para- 
metrically iso-spectral, that is da{£)/dtj — for all tj G M. The iso-spectrality condition gives rise 
to a hierarchy of commuting to each other nonlinear bi-Hamiltonian dynamical systems on the 
functional manifold M in the general form 

(2.3) ^{u{t),v{t)Y = -i9gradH,[u,v] := K,[uit),v{t)], 

where Kj : M ^ T{M) and Hj E 'D{M),j E Z+, are, respectively, vector fields and conservation 
laws on the manifold A/, which were before described in [Sj IH [7] , 

d 
d 



(2.4) d 



is a Poisson structure on the manifold M and, by definition, 
^2 5) I (u{x,ti,t2tz..:) 

for t E M"^ 



It is well known [151 [HI E [S] that the Casimir invariants, determining conservation laws for 
dynamical systems (j2.3p . are generated by the suitably normalized monodromy matrix S{x; A) G 
End of the linear problem 



(2.7) —F{y,x;X)^e{y;X)F{y,x;X), F{y,x; X)\y=, ^ I, 



(2.6) S{x; A) = k{X)S{x; A) - ^trS'(a;; A), 

where F(y, x; A) G End <C'^ is the matrix solution to the Cauchy problems 
d_ 
dy' 

for all A G C and x,y E M., where I EEnd is the identity matrix, S(x; A) :— F{x + 27r, x] A) is 
the usual monodromy matrix for the equation (j2.7p . Here the parameter k{X) G C is invariant 
with respect to fiows (12.31) and is chosen in such a way that the asymptotic condition 

(2.8) S{x-\)Eg- 

as A — > cxD holds for all a; G M. Here C where Q :— Q+ ® Q- is the natural splitting 

into two affine subalgebras of positive and negative A-expansions of the centrally extended [THl [2S] 
affine current s[(2)-algebra Q := Q ®€- : 

(2.9) Q:^{a= ^ a^^^ ® X' : a^^'> G C°° {R/2nZ; s((2; C)) }. 
The latter is endowed with the Lie commutator 

(2.10) [(ai,ci), (a2,C2)] := ([ai, 02], (ai, da2/(ia;)), 
where the scalar product is defined as 

(2.11) (ai, 02) := resA=oo / tr{aia2)dx 

Jo 

for any two elements 01,02 G G with "res" and "tr" being the usual residue and trace maps, 
respectively. As the spectrum a{£) C C of the problem (|2.2p is supposed to be parametrically 
independent, fiows (|2.3I) are naturally associated with evolution equations 

(2.12) dS/dtj = [{X^+'^S)+,S] 

for all j E M, which are generated by the set I{Q*) of Casimir invariants of the coadjoint action 
of the current algebra 5 on a given element £{x] A) G Q*_ = G+ contained in the space of smooth 
functionals T>{Q). In particular, a functional 7(A) G I{G) if and only if 

(2.13) \S(x; A), £(x; A)] -I- -^S{x; A) = 0, 

dx 
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where the gradient S{x; A) := grad7(A)(^) G is defined with respect to the scalar product (|2.11l) 
by means of the variation 

(2.14) 57(A) (grad7(A)(£),J£). 

To construct the solution to matrix equation (12.131) . we find preliminary a partial solution 
F{y,x;X) G End C^, a;, y G R, to equation (|2.7|) satisfying the asymptotic Cauchy data 

(2.15) F{y,x;X)\y=,=I + 0{l/X) 
as A cxD. It is easy to check that 



(2.16) Fiy,x;X) = 



ei(y,x;A) ~^%^e2(y, x; A) 
^^ei{y,x;X) e2{y,x;X) 



is an exact functional solution to (|2.7p satisfying condition (I2.15|) . where we have defined 



y 



(2.17) ei(?/,x; A) := exp{^[u(a:) - it(?/)] + A / a dv{s)}, 



y 



e2iy,x;X) ■.= exp{^[u{y) - u{x)] - ^ J l3 ds}, 

with the vector-functions G C°°(R/27rZ; M) satisfying the following determining functional 
relationships: 

(2.18) /3 = -(u2-2«,+C/3)'^', 



as ^ := 1/A — > and existing when the condition (p{x, t) :— ^Ju\ — 2vx ^ on the manifold M at 
t = G M'^. 

The fundamental matrix F{y, x; A) G End can be represented for all y G M in the form 

(2.19) Fiy,x;X) = F{y,x;X)F-\x,x;X). 
Consequently, if one sets y = a; + 27r in this formula and defines the expression 

(2.20) k{X) := X-^[ei{x + 2TT,x;X)~e2ix + 2TT,x;X)]-^, 

it follows from p.6p . (|2.16p and ()2.19p that the exact functional matrix representation 



la(xi\)+l3{x:X)] 

2X[a(x:X)-l3{x;X)] \^[a(x;\)-l3(x;\)] 

1 _ [P(x;X)+aix-X)] 

[a{x;X)~f3{x;X)] 2X[l3{x;X)-a{x;X)] 



(2.21) S{x;X) 
satisfies the necessary condition (|2.8I) as A -> 00. 

Remark 2.1. The invariance of the expression (|2.20p with respect to the generating vector field 
(|2.3p on the manifold M derives from the representation (|2.19l) . the equations (|2.13p and 

d X^ 

(2.22) -F{y, xo; m) = J^S{x; A)F(y, xq; m), 

which follows naturally from the determining matrix flows ()2.12p upon applying the translation 
y ^ y + 27r. 

The matrix expression (I2.2ip gives rise to the following important functional relationships: 
1 - A(§ii - §22) . -2A^.§i2 3 

2S21 1-A(sii-S22) 

which allow to introduce in a natural way the vertex operator vector fields 
(2,24) .Yj.cxp(±D,), Z,, 2 

acting on an arbitrary smooth function 77 G C°°(M^+;M) by means of the shifting mappings: 
(2.25) 



T](x,ti,t2,...,tj,...) := T]^{x,t;X) 



= r^{x,h ± 1/A,i2 ±/(2A2),t3 ± l/(3A3)...,t, ± l/(jAJ),...) 
as A — )■ 00. Namely, we following proposition holds. 
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Proposition 2.2. The functional vertex operator expressions 

(2.26) a{x,t;X) = X^a{x,t) ^ {x,t;X), 

/3(x,i;A) = X+p{x,t))^p+{x,t-\) 

solve the functional equations V2.18\) . that is 

a- = + {ul -2v^+^a-Y/^, 

(2.27) p+ = - {ul - 2v, + 
where t G and ^ = 1/A 0. 

Proof. To state this proposition it is enough to show that the fohowing relationships hold: 



(2.28) 



d 
d_ 



1 - A(sii - S22) 



2521 



1 - A(sii - S22) 



A=l/? 



A=l/C 



d 

di 

d_ 
'di 



1 - A(sii - S22) 



2S21 
'SA^Sia 



for any parameter ^ — 0, where by definition 

d d 



(2.29) 



dt ■ d^ 



1 - A(sii - S22) 



A=l/« 



A=l/C 



dt 



is a generating evolution vector field. Before doing this wc find the evolution equation 
(2.30) 



lsix■,^,)^[X^^Six■,^,),S{x;X)] 



on the matrix S{x; fi) as /i, A — > 00, which entails the following differential relationships: 



(2.31) 



dsii/dt = A^(s2ic?si2/dA - si2rfs2i/rfA), 
ds22/dt = X^{si2ds2i/dX — S2idsi2/dX), 
dS22/dt = A3[Si2^(sii - S22) - (Sii - §22)^), 



dSiildt = A3[S2ii^(S22 - hi) - (S22 - hi)^)- 

Using these relationships p.3ip . one can easily obtain by means of simple, but rather cumbersome 
calculations, the needed relationships (|2.28p . As their direct consequences the vertex operator 
representations (|2.26p for the vector functions a, /3 G C(M^+ ; R) hold. □ 

Now we take into account that, owing to the determining functional representations (|2.18p . that 
the limits°° 



(2.32) 



lim a (x,t;X) = Ux{x,t) + ip{x,t), 

A—)- 00 



lim /3+(a;,i; A) = Ux{x,t) - ip{x,t), (p{x,t) := ^u\(x,t) - 2vxix,t), 



exist on the manifold M. Moreover, having iterated the functional relationships (|2.18p . one can 
find that 

X.a = a =Ux + ip + E,{ 1 ) + 

LP 

. "L + '^Uxx^Px - U-ix^ ^ ifxxV + 5(/?g ^ ^ 3 

2 p^ ipfl ^ 



(2.33) 



Xtfi - /3 



ip ip 



^ ul^ - 2Uxx(Px + Uzx^ ^xx^ + 5t/?g , 
2 Lp^ Lp^ ' 



which immediately yield the higher Riemann type commuting nonlinear Lax integrable dispersive 
dynamical systems on the functional manifold M. For instance, making use of the relationships 



(2.34) 



lim [a"(x,i;A) ± ;3+(x, i; A)]/2 

A— >-oo 



Ux{x,t), 

(p{x,t) 
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one easily obtains that 



(2.35) , 



and so on, where ip — \Jv?^ — 2vx and we took into account that the following asymptotic expan- 
sions hold 



(2.36) 



X^Kx,t;X) = Ux -ip + ^{ux,ti ^ V'ti)+ 

+ ^(Ux,tuti - VUM + UxM - Vts) + O(C^) 

as C = 1/A 0. 

It is worth here to mention that the scheme devised above for finding the corresponding vertex 
operator representations for the Riemann type equation (|2.ip can be similarly generalized for 
treating others equations of the infinite hierarchy (jl.ip when N > 3, having taking into account 
the existence of their suitable Lax type representations found before in recent works [Ml [51 [3] • 

3. Concluding remarks 

The vertex operator functional representations of the solution to the Riemann type hydrody- 
namical equation (|2.ip in the form (|2.27p is crucially based on the representations (|2.23p and 
evolution equations (j2.28p . which provide a very straightforward and transparent explanation of 
many of "miraculous" vertex operator calculations presented before both in [T31 [2] and in . It 
should be noted that the effectiveness of our approach to studying the vertex operator represen- 
tation of the Riemann type hierarchy owes much to the important exact representation (j2.2ip for 
the corresponding monodromy matrix, whose properties are described by means of applying the 
standard [15, 18][Tni[10] Lie-algebraic techniques. As an indication of possible future research, it 
should also be mentioned that it would be interesting to generalize the vertex operator approach 
devised in this work to other linear spectral problems such as those related to dynamical systems 
with a parametrical spectral (TI] [T71 TS' dependence, spatially two-dimensional ^lOj , Pavlov's and 
heavenly [12. dynamical systems. 
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